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Abstract 

We outline a partial-fractions decomposition method for determining the one-particle spectral 
function and single-particle density of states of a correlated electronic system on a finite lattice 
in the non self-consistent T-matrix approximation to arbitrary numerical accuracy, and demon- 
strate the application of these ideas to the attractive Hubbard model. We then demonstrate 
the effectiveness of a finite-size scaling ansatz which allows for the extraction of quantities 
of interest in the thermodynamic limit from this method. In this approximation, in one or 
two dimensions, for any finite lattice or in the thermodynamic limit, a pseudogap is present 
and its energy diverges as T c is approached from above; this is an unphysical manifestation 
of using an approximation that predicts a spurious phase transition in one or two dimensions. 
However, in three dimensions one expects the transition predicted by this approximation to 
represent a true continuous phase transition, and in the thermodynamic limit any pseudogap 
predicted by this formulation will remain finite. We have applied our method to the attractive 
Hubbard model on a three-dimensional simple cubic lattice, and find, similar to previous work, 
that for intermediate coupling a prominent pseudogap is found in the single-particle density of 
states, and this gap persists over a large temperature range. In addition, we also show that for 
weak coupling (an on-site Hubbard energy equal to one quarter the bandwidth) a pseudogap 
is also present. The pseudogap energy at the transition temperature is almost a factor of 
three larger than the T=0 BCS gap for intermediate coupling, whereas for weak coupling the 
pseudogap and T=0 BCS gap energies are essentially equal. These results show that a pseudo- 
gap due to superconducting fluctuations occurs in three dimensions even in the weak-coupling limit. 

Keywords: Correlated electrons - theory, Hubbard model, spectral functions 
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I. INTRODUCTION 



The theory of Bardeen, Cooper, and Schrieffer (BCS)^ has been very successful in ex- 
plaining a wide variety of phenomena associated with superconductivity. This is all the more 
remarkable since it is based on the independent quasi-particle approximation. Since very 
early on corrections to this simple picture have been investigated to explore what effect, 
if any, the breakdown of the quasi-particle approximation would have on various super- 
conducting properties. One body of work in this category is the investigation of collective 
excitations, beginning with Anderson,— Bogoliubov et a/.—, Rickayzer4 and others.— Much of 
this work was directed towards a proper gauge- invariant description of the Meissner effect. 

A separate line of work started with Migdal^ and Eliashberg — and investigated how the 
coupling of electrons to dynamical phonons would affect the quasiparticle properties of the 
electrons, including the transition temperature. 

While this and other work 8 focussed on properties in the superconducting state, a dif- 
ferent line of enquiry concerned pairing fluctuations above T c , as these could possibly cause 
deviations from the mean field picture provided by Fermi Liquid/BCS theory. Within BCS 
theory superconductivity arises through the formation of an order parameter, which has both 
an amplitude and a phase.— Fluctuations can impact both the amplitude and the phase. At 
the present time phase fluctuations have been suggested to be relevant to the high-T c prob- 
lem. In particular, a school of thought 10 ' 11 ^ 2 ! 13 suggests that electrons have a pair amplitude 
already above T c , i.e. they exist as 'pre-formed' pairs. The superconducting transition is 
then simply the temperature at which the pairs acquire a common phase.— The extreme case 
of this scenario is where two electrons already form a well-defined boson at the critical (i.e. 
phase-forming) temperature; this temperature is then the Bose condensation temperature 
for the system. A comprehensive review of such ideas has been provided by Randeria.— 

These so-called pairing fluctuations would generally lead to a non-zero amplitude for 
pairing above T c , with long-ranged superconducting order destroyed by fluctuations. One 
school of thought (to which we subscribe) suggests that some of the anomalous normal state 
properties observed in the high T c systems, in particular the pseudogap, are due to pairing 
fluctuations, and a goal in this paper is the examination of such fluctuations in different 
dimensionalities using one particular controlled approximation scheme. 

Specifically, one technical framework with which one can hope to understand this physics 



3 



is the so-called T-matrix approach. Within such a scheme, a pairing susceptibility can be 
defined;— in the usual case, 17 this pairing susceptibility diverges at T c , the latter defined 
by BCS theory. Finding this temperature constitutes an alternate means of defining the 
transition temperature, and is generally known as the Thouless criterion.— However, the 
existence of a diverging susceptibility (as one approaches the BCS critical temperature from 
above) should impact the single electron properties that determine the critical temperature 
in the first place. The T-matrix approximation accomplishes this through a renormalization 
of the potential an electron experiences. However, a limitation of this approach is that it 
can only be justified as a low density expansion of the two-particle Green's function . 19 i 20 

Several many-body frameworks within which the T-matrix approximation can be cast 
have been discussed by Baym and Kadanoff and coworker s 21 i 22 i 23 many years ago. The 
particular choice of framework unfortunately influences the results one obtains when cal- 
culations for a particular model system are actually performed. For example, in Ref.,— a 
particular prescription for the single electron self-energy is derived through an equation-of- 
motion method. The resulting expression requires the use of two 'levels' of single electron 
Green functions, one bare, the other renormalized; that is, the latter contains a self-energy, 
and, consequently, such calculations must proceed self-consistently. Nonetheless, already in 
this work 2 ^ (see footnote 13 and the discussion at the end of section 3) it was recognized 
that certain features of this approximation were not acceptable. 

Several years later the issue was debated some more , 24 > 25 i 26 until what one could refer to 
as the 'Patton fix' emerged as the allegedly correct way to handle the difficulties with the 
Kadanoff-Martin prescription. In this adjustment, one of the self-consistent single electron 
Green functions is simply replaced with its non-self-consistent bare counterpart. This cor- 
rects some problems, but the adjustment is completely ad hoc, and has been called "non 
derivable" Such an approach has been championed recently by Levin and coworkers,— 
and has found some agreement with experiments on the high T c cuprates. 

Earlier, Baym^ 3 - had formulated a theory of single electron and two particle properties 
based on functional derivatives of a free energy functional, a procedure which guaranteed 
that certain conservation laws would be satisfied. This formulation is known as a conserving 
approximation, and its extension to lattice electrons is known as the FLuctuation Exchange 
approximation (FLEX) 29 . When only the particle-particle channel is retained, the FLEX 
becomes a T-matrix theory with fully self-consistent single electron Green functions. Such 
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a formulation would appear to be the most accurate of the T-matrix approximations, since 
it includes the most number of diagrams. However, this claim is, in our opinion, unfounded. 
The possibility that vertex approximations (omitted in all T-matrix theories discussed here) 
partially cancel some of the self-energy contributions retained in the more fully self-consistent 
T-matrix approximations appears likely— That is, including more channels and making 
everything self-consistent is not the panacea it might seem to be — this point has been 
emphasized in the alternative formulation of this problem by Vilk and Tremblay.— 

It is probably fair to say that 'modern' interest, and a viable theoretical approach, on the 
influence of pair fluctuations on the superconducting transition was initiated by Leggett.— 
This was later followed by an important paper by Nozieres and Schmitt-Rink,- 33 « wherein they 
solved for the zero temperature properties of an electron gas with attractive interactions, for 
all coupling strengths. They showed that the BCS approximation gave the correct limiting 
behaviour in the strong coupling limit, where, in the lattice model at least, it was clear that 
the model was best described by tightly bound pairs of electrons. In this limit, it became 
clear that the superconducting transition was properly described by a Bose condensation of 
electron pairs. An important technical point raised in these two papers was that both the 
BCS gap equation and number equation were required, and had to be solved self-consistently. 

This work was further extended to two dimensions in the context of the high T c cuprates, 
by Schmitt-Rink, Varma, and Ruckenstein, in Ref.— They made the point that, for a specific 
(and fixed) electron density, as the temperature was lowered the chemical potential would 
self-consistently adjust so as to avoid the Thouless instability. (Note the necessity of adhering 
to the lesson from Ref.,— that being that one must solve both for the normal state properties 
and the electronic density simultaneously.) In fact, all electron densities would result in a 
situation at T = in which the Fermi surface no longer exists; the electrons are paired in 
bose-like bound states, and the chemical potential is simply half of the single pair bound state 
energy. This result violates Luttinger's theorem.^ Improvements were made by Serene,^ 
Tokumitu et al£L, and later, for lattice fermions, by three of the present authors,-^ but the 
physics obtained in Ref.-*^ remains. In fact, for lattice fermions the results are more physical, 
as the Thouless curve becomes a locus of points in the T — \x plane on which the electron 
density is identically equal to unity 38 (as opposed to a diverging quantity). Again, the 
simultaneous solution of the superconducting transition and the equation for the electronic 
density is required. 
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Since then a significant number of calculations have appeared on the T-matrix prob- 
lems; see, for example, Refs . 39 i 40 i 41 i 42 i 43 i 44 i 45 i 46 i 47 In addition, arguments against such for- 
mulations and in favour of alternate frameworks requiring vertex corrections, have also 
appeared . 31 i 48 i 49 i 50 This body of work represents a significant advance in the study of the 
effect of fluctuations at temperatures above the superconducting transition. However, much 
of the 'mathematical contortions' that are invoked can be traced to the difficulty of studying 
two-dimensional systems in which the Mermin- Wagner theorem precludes the existence of 
true long-ranged order (in the systems studied here).— Self consistency, and other formula- 
tions, aim to recover this result from any Thouless criterion approach (by which we imply 
theory attempting to identify T c through a diverging pair susceptibility). 

A simple remedy to avoid the Mermin- Wagner theorem is to study a three dimensional 
system. Therefore, in this paper we focus on a purely electronic Hamiltonian that yields a 
nonzero superconducting transition temperature, the attractive Hubbard model for a sim- 
ple cubic lattice. Further, in three dimensions we can avoid the necessity of using a self- 
consistent theory, especially for weak coupling, and this is a simplification which will allow 
for us to solve for the self energy of a finite lattice, using the method outlined in the following 
sections of this paper, to essentially any numerical accuracy that we require. 

In the above-described work, the appearance of a pseudogap, an experimentally observed 
phenomenon in which the single-particle density of states is suppressed in the form of an 
incompletely formed superconducting gap,— is predicted. The literature on this question 
is now substantial — for a review of recent ideas see, for example, Ref.-2£ A recurring 
theme in some of this work (see, e.g. Ref— -) is that the pseudogap found in the high T c 
cuprates can be associated with the two dimensional nature of the CuO planes. In terms of 
fluctuations, and in its simplest form, this is simply the statement that one expects enhanced 
fluctuations, towards the low-temperature phase, the lower the effective dimensionality of 
the system. The analysis that we present in this paper suggests that two dimensions may, 
in fact, not be necessary, because a pseudogap is found in the three dimensional attractive 
Hubbard model at weak coupling. For example, even at a temperature of a few per cent 
above the superconducting transition, in the three-dimensional attractive Hubbard model 
at half filling, clear evidence of a single-particle density of states pseudogap is found, and 
the magnitude of the gap is comparable to the BCS gap that would be found at T = 0. 

We note that the number of papers on the 3d pseudogap is very small, compared to that 
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for the 2d problem, so comparisons with other work are difficult. Instead, here we shall 
simply address the question, should one expect a pseudogap in the 3d system with purely 
electronic attractive interactions? 

II. FORMULATION 

A. Model 

We consider lattice fermions and the attractive Hubbard model; we restrict our attention 
to a kinetic energy that includes nearest neighbour hopping only. Working in the grand 
canonical ensemble thus leads us to consider the so-called Grand Hamiltonian /C defined by 

K. = -t ( C Lx C m',a + C ] mlu C m ^) -fl n ™ ~ \ U \ E n mt n ml • (1) 

(m,m') m,cr m 

In this equation t is the single particle hopping matrix element for nearest neighbour sites, 
H is the chemical potential, \U\ is the on-site attractive potential, m and m' label the lattice 
sites, while (m, m') denotes nearest neighbour sites (which are counted only once in the 
kinetic energy sum), and o is the electron spin. The electron creation (annihilation) and 
number operators with spin a at lattice site m are given by cj^ (c mcr ) and n ma , respectively. 
We restrict ourselves to one-dimensional chains, and three-dimensional simple cubic lattices, 
and employ periodic boundary conditions. 

B. One and Two Particle Propagators 

Firstly, we define the required thermal Green's functions. The single electron Green's 
function^ monitors the history of an electron injected at (imaginary) time t[ and lattice 
site m[, and removed at time t% and site mi, and is defined by 

G T ( mi , Turnip) = (-)(T T [c miT (ri)i, T (rO]) = G T (11') = G(l - l') . (2) 

We use standard abbreviating notation, 1 = (m 1; ti), and note that the last result of Eq. (J2J) 
follows from the translational periodicity, the time invariance, and spin isotropy of our lattice 
and Hamiltonian. Similarly, the two-electron propagator is given by 

^Tl(mi,n;m2,r 2 ;m' 1 ,r{;m2,r2) = 0(12; 1'2') 

= (TV [c miT (ri)c mai (r 2 )ci 1 , <L (7^)ci lit (^)] ) . (3) 
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(In this paper we will need to only consider the function, viz. the singlet channel, and thus 
spin indices of Q will not be shown.) This latter function's transform in momentum/energy 
space can be simplified by incorporating translational periodicity and temporal invariance, 
and we introduce the notation Q = (Q, %vq) and Ql = Q • mi + %vqT\. Then we have 

0(12; 1Y) = — ^ ^ e 4(a+2)-d'+2')) e ^d-2) e W(i'-2') g(Q ;qq >) (4) 

(PN) Q qql 

where N is the number of lattice sites, Q is the centre-of-mass generalized momentum of 
either the outgoing or incoming pair, and q (q') is the relative momentum of the outgoing 
(incoming) pair of electrons, respectively All wavevector summations span the first Brillouin 
zone, and Matsubara sums go over all integers. 

From this function we extract the quantity that is our principal concern, viz. the pair 
propagator QP air which can be included in this set of equations in the following two ways: 

G pair (Q) = 7^-2 £ G(Q;qq') (5) 

lP iV ) qq' 

£(n + ; = j-^j £ e^-^g^iQ) . (6) 

An examination of Eqs. (JUEJ) shows that the pair propagator is equal to the probability of 
adding a pair of (opposite spin) electrons with one particular centre-of-mass momentum Q 
to one point in space, and then removing them as a pair at some other point in space (the 
centre-of-mass momentum Q is conserved at all collisions that occur between these points 
in time, and one must sum over all possible relative momenta, q and q'). 

The pair propagator's connection to superconductivity can be made apparent by intro- 
ducing the pair annihilation operator via 



A Q = c § +qT c a- qi (7) 

and then noting, by Eq. (j5J), that the pair propagator can be expressed as 

G pair (Q) = f o e^(T T [A Q (T)A^(0)]> air . (8) 

This equation most clearly allows one to understand the association of the pair propagator 
with the so-called pair susceptibility.^ 

Finally, it is instructive to extrapolate these relations to real times — one sees that one 
other way to interpret the behaviour of the pair propagator is that its analytic continuation 
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from imaginary frequencies to the complex plane allows one to evaluate the retarded pair 
propagator, in real time. As the discussion of Ambegaokar— makes apparent, the poles of 
the pair propagator correspond to the relaxation times of a pair described by the retarded 
Green's function (also see Ref.—). Most of the poles correspond to continuum states com- 
prised simply of two single (quasi-) electron states. Two poles are exceptions; initially these 
occur on the real axis. However, at sufficiently low temperature these poles form conjugate 
pairs on the imaginary axis. One of these poles signals an exponentially unstable (viz. grow- 
ing) paired state implying that the normal state has now become unstable with respect to 
some new state whose stability is presumably related to these paired electrons. It is clear 
that, using this approach, one cannot proceed to temperatures below the critical tempera- 
ture T c . On the other hand, studying the pair propagator above this critical temperature is 
an ideal way to study pair fluctuations, as well as other properties of the fully interacting 
normal state. This formalism for identifying T c from above has come to be known as the 
Thouless criterion.— 



C. Non Self-Consistent T-Matrix Approximation 

The above discussion motivates the analysis of the pair propagator in the normal phase. 
In order to find this function for the AHM we need to invoke a number of approximations. 
To be specific, we employ the ladder approximation to the Bethe-Salpeter equation for the 
reduced vertex. The derivations associated with this approximation, which leads to the T- 
matrix approximation for the single-particle propagator, are well known. It has been shown 
that this approximation is justified (see, e.g., Ref.—) in the limit of vanishing electron (or 
hole) densities. 

The physical content of this approximation can be made most apparent if we begin with 
the equation for the 2-particle Green's function expressed in terms of the effective 2-particle 
interaction, which is given by (again, suppressing spin indices and sums over these indices) 

G(Q- q, q') = (3N6(q - q')G°(^ + g)G°(| - q) 
- + g)G°(| - q)V(Q; q, g')G°(| + g')G°(f - q') (9) 

where G° is the noninteracting single-electron Green's function. The ladder sum for the 
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effective interaction of the AHM, in a non self-consistent (NSC) calculation, is given by 

r(Q; q, q>) = - \U\ + \U\ ]TG°(| + fc)G°(| - k)T(Q; k, q') . (10) 

k 

This demonstrates that in the ladder approximation one includes scatterings of all orders 
in the particle-particle channel. Further, due to the independence of V on internal relative 
momenta, the solution of Eq. (flT)|) is straightforward, viz. 

where the pair susceptibility of the noninteracting system is equal to 

X°(Q) = ^EG°(| + g)G (|-g). (12) 

Finally, using Eqs. ()5p9Jllj) we can derive the pair propagator, or pair susceptibility x(Q), 
in the NSC ladder approximation: 

= *«> = ; _g (0) (13) 

The solution of the Thouless criterion for the AHM treated in the NSC ladder approxima- 
tion is now clear. The lifetime of a newly introduced (singlet) pair of electrons with crystal 
momentum Q will first diverge (on cooling) when Imx°{Q) = and 1 — \U\x°{Q) = 0. 
The analytical properties of the pair susceptibility dictate that this must occur for zero 
Matsubara frequency, namely the Thouless criterion reduces to 

X°(Q,0) = . (14) 

Hereafter, the locus of points in the T — \i plane which satisfies Eq. (|14jl shall be referred 
to as the Thouless criterion curve. For the AHM it can be show n i 55 that on cooling the 
transition defined by Eq. ()14|) first occurs for Q = 0, and from now on we define T c in this 
manner. 

As detailed in the introduction, we are interested in gaining an understanding of the effects 
of 'feedback', 56 that is when the pair fluctuations that are included in the pair propagator 
are introduced into the description of the single-electron properties. To this end, we evaluate 
the proper self-energy in the (same) NSC T-matrix approximation, resulting in 

E(*)= -JW^ X (Q)G°(Q-k). (15) 
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(Note that throughout this report we ignore the Hartree term, which can be simply absorbed 
into the chemical potential.) From this quantity the fully interacting single-particle Green's 
function is found via Dyson's equation 

G{k) = - m ■ (16) 

The important quantity in this paper, the single-particle density of states (DOS), is related 
to the one-electron spectral function, given by 

A(k,u) = - - ImG(k,cu + iO + ), (17) 

which leads to the DOS given by 

A/» = ££A(k,u;). (18) 
D. Partial Fractions Decomposition 

The above equations are well known. Here we present an approach that allows for us 
to solve for all of the spectral functions of finite lattices in the above-described NSC T- 
matrix approximation to arbitrary numerical accuracy — in one dimension we can solve 
for chains of lengths up to 128 (which has 65 inequivalent k points in the first Brillouin 
zone), while in three dimensions we have solved lattices up to a size of twelve cubed (which 
has 84 inequivalent k points in the first Brillouin zone) — and to then extract information 
about the density of states in the thermodynamic limit. To complete the decomposition for 
a twelve cubed lattice at one temperature and chemical potential can take up to ten days 
on a high-speed 64-bit (which is essential) alpha processor; so, it is not a simple matter to 
obtain large sets of data using this approach. 

To solve for the one-particle self energy we first note that the non-interacting pair sus- 
ceptibility, defined in Eq. p2jl . can be written in terms of a sum over wave vectors 

x°(Q.^.,) = ^E / l kl+ f Q ; kl ' 1 - < 19 > 

where /[£] represents the Fermi-Dirac distribution function at a temperature T and chemical 
potential £k = £k — fi, and is the noninteracting band electron dispersion (e.g., = 
— 2tcos(k) in one dimension). The analytic continuation of this function to the complex 
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plane, denoted by x ~ > Xi is achieved trivially by the substitution iv n — > z. Thus, from 
Eq. (fE?j) we can write the pair propagator as 

If one simply expands the wave vector sums, it is seen that this function admits a partial 
fraction decomposition. Incorporating the requisite analytical properties of this 2-particle 
propagator (the chosen signs follow from the spectral representation of the 2-particle GF) 
we choose to write this expansion as 



X(Q, z) = - K Q J 1} Q - 1 Q J 2) Q - • • • (21) 



where the number of poles for each Q is denoted by sq, the energies Eq are real, and the 
residues Rq are strictly positive, for £ = 1, ...,sq. Substituting this into Eq. (JT5j) and 
completing the Matsubara frequency sum analytically allows us to restate the expression for 
the self energy in the convenient form 

where b[E] is the Bose-Einstein distribution for temperature T with the chemical potential 
set equal to zero. 

The application of the above results requires the determination of the poles and residues 
of the partial fraction expansion of the pair propagator. This is a complicated task, and we 
have completed our calculations using a computer algebra system allowing for a symbolic 
calculation that produces numerical results for the poles and residues. 57 Also, note that 
the result of this partial-fractions formulation is a DOS that is a series of delta functions - 
we have broadened each of the single-particle spectral peaks with Lorentzians with a fixed 
broadening of 3% of the bandwidth, namely 0.12£ in one dimension, and 0.36t in three 
dimensions. 
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III. EXAMINATION OF THE ID PSEUDOGAP - TEST CASE 



As an example of the application of this method, we have examined one dimensional (Id) 
chains. We stress that this is not, and is not intended to be, a solution of the Id AHM, 
which cannot be accomplished through the use of the T-matrix approximation. Instead we 
use these results as a means of demonstrating how our results can be applied, via a finite- 
size scaling extrapolation to the thermodynamic limit, to estimate the bulk value of the 
(pseudo)gap. 

In Fig.^we show the DOS of a one-dimensional chain with periodic boundary conditions 
of varying lengths (16 x 1 up to 128 x 1) with \U\/t = —2 in the non self-consistent T- 
matrix approximation for the AHM. We have taken this data for a fixed (non-interacting) 
electronic density per lattice site of (n) = 0.5, corresponding to 1/4 filling. The temperature 
is fixed to be just above the BCS transition temperature, and to account for the fact that the 
BCS transition temperature depends on the length of the system, we have used a varying 
temperature of 1.001 T C (L), thus corresponding to a fixed ratio of T/T C (L) for all system 
sizes. 

It is clear from Fig. [T] that each DOS curve contains sharply peaked structures, and it is 
natural to associate the location of these peaks with something analogous to the supercon- 
ducting gap. Of course, these data are evaluated in the normal phase, and the DOS does not 
ever go strictly to zero (near ui/t = 0) in the normal phase, so the locations of the peaks are 
not true superconducting gaps, but can be thought of as something akin to (single-particle 
DOS) (pseudo)gaps. 

These data display strong finite size effects, and the obvious question is, what does the 
DOS look like in the thermodynamic limit? We now demonstrate how one can answer this 
question, and thus use data such as that shown in Fig. ^ to correctly extrapolate to the 
thermodynamic limit. This analysis will lead us to a controlled manner of both verifying 
the existence of, and quantitatively estimating, the pseudogap in three dimensions. 

The finite-size scaling corrections to the wave vector sums that are encountered in this 
problem go as 

i? = s/> + (i?-s L iq ) ~ I L ^ + °^»- (23) 
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FIG. 1: The density of states (DOS) of the Id attractive Hubbard model evaluated in the non 
self-consistent T-matrix approximation for 16x1, 32x1, 64x1, and 128x1 chains with periodic 
boundary conditions, for \U\/t = —2 and a fixed (noninteracting) electronic density per lattice 
site of (n) = 0.5, corresponding to 1/4 filling; the temperature used for each curve is 1.001 T C {L). 
Each delta-function peak of the spectral function has been broadened as a Lorentzian using a fixed 
broadening energy of 3 % of the bandwidth — through experimentation we have found that such 
a broadening produces smooth density of states curves without washing out important structure. 



Thus, noting that the peak locations, which from now on are provocatively labelled as A, 
from Fig. ^ can be associated with a strong divergence of one particular contribution (see 
below) to the self energy, and that the wave vector sums in Eq. ()15j) should have the same 
finite-size scaling corrections as given in Eq. ([23)1 . we expect that the peak locations of, e.g., 
Fig.^ should approach their values in the thermodynamic limit as ~ 1/L 2 . This is precisely 
the form that we find, as is shown in Fig. 121 

We now focus on the values of A obtained from using the above-shown finite scaling 
analysis extrapolated to the thermodynamic limit. To be specific, we now show that the 
divergence of this (pseudo)gap as T approaches T c from above is precisely what should be 
found in the thermodynamic limit. This successful extraction of A(L — > oo) helps to support 
our contention that such a finite-size scaling analysis is a reliable procedure for extracting 
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FIG. 2: The locations of the peaks, labelled as A and expressed in units of t, found in the Id DOS 
(e.g., see Fig. ^), plotted as a function of one over the the system's length (L) squared, for a variety 
of temperatures (shown in the legend) just above the superconducting transition temperature. The 
same parameters as in Fig. n are used. The solid lines are least squares fits to the data for each of 
these temperatures. 

the bulk value of this energy. 

To this end, we first require the behaviour of the self energy as the Thouless criterion is 
approached from above. It is instructive to return to the definition of the self energy (prior 
to making a partial fractions expansion) given in Eq. (fTK|) . All quantities in Eq. (|TH|) are 
known exactly (although X°(q>2z/ n ) requires a momentum sum - see Eq. ()12j)); nonetheless, 
it is not in a convenient form from which the DOS in the L — > oo limit can be extracted. 
That is, while one can successfully approximate the wave vector integration over the first 
Brillouin zone with discrete momentum sums (with a dense mesh of k points) for higher 
temperatures, such a method is doomed to fail for temperatures/chemical potentials near 
the Thouless curve. This is because the denominator in the kernel of Eq. Ijl5j) (i.e., the 
denominator of Eq. ([13))) approaches zero (for q = and iv n = 0). To be specific, using a 
discrete momentum sum over the wavevectors of a finite lattice will result in a self energy 
that diverges as T — > T c in any dimension, for any finite lattice. However, phase space 
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arguments following the Mermin- Wagner theorem^ show that the self energy in the form 
of Eq. (fT3j) should not diverge in three dimensions - see below. Thus, we need to examine 
Eq. (JT5J) more carefully to properly understand the distinction between one (and two) and 
three dimensions, and to properly extract the temperature dependence of the (pseudo)gap 
as T c is approached from above. 

Since the Thouless criterion arises from the %v n = term in the Matsubara frequency, it 
is prudent to treat the n = part of the self energy sum separately. Then, to correctly deal 
with the q — > contribution, we have separated the %v n = term into two contributions - 
the one from q ^ is given by 

E nzq (k,iu; m ) = -\U\ 2 ^J^ X (^0)G°(-k + q,-iu m ). (24) 

and the one from q near zero is done by an integration (which is correct, and as we show 

necessary, if we completed the entire calculation in the thermodynamic limit), 

1 rQo d^n 1 

£ 29 (k^ m ) = -|f/| 2 -/ 7^x(q,0)— -, (25) 

P J-qo (2n) d -iuJ m ~ (e-k+q - AO 

where qo represents an upper cutoff for the part of the first Brillouin zone in which we will 

integrate over the integrand analytically, and d represents the dimensionality of the space. 

The prime in the summation of Eq. (}2*4*j) indicates that this same small q region is omitted 

from Tj nzq . Thus, in all expressions except the last (Eq. (|25|l ) we will simply perform q-sums 

over finite lattices (which amounts to the trapezoidal rule for numerical integration), while 

in the latter equation it is important to integrate analytically over the integrand. 

To this end we require an explicit expression for x (q, 0), which is difficult to obtain. 

Instead, we approximate it with an expression valid for small q: 

X °(q,0)^a 2 -6V, (26) 

where q = |q| and a 2 and b 2 are obtained by fits to x°(q, 0). In one dimension this is a simple 
matter of fitting the susceptibility at q = and at some small q — qo, and an example case 
is shown in Fig. |21 clearly, this is an excellent approximation for small q. Then, expanding 
the band dispersion e^k+q for small q, the integral in Eq. (|23|) can be written as 



(27) 
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FIG. 3: Real part of the non-interacting susceptibility vs. q in one dimension. Red squares are for 
a 16 site lattice, while the green curve is based on the fit as described in the text. The blue curve 
is the result in the thermodynamic limit. Clearly, the black curve reproduces very accurately the 
thermodynamic limit result at low q, even though the green curve is obtained on the basis of a fit 
to the two lowest points (q = and q = n/8). 

We have introduced A = 1 — \ U\a 2 and B = \U\b 2 , and note that both are positive numbers. 
Further, and most importantly, since x°(0, 0) ~ 1/T at low temperatures, we obtain that A 
is approaching zero linearly in T from above asT-> T c . That is, A ~ |a|(T — T c ) for alpha 
being some constant. (This clarifies the "mean-field" label that can be applied to the non 
self-consistent T-matrix approximation.) 

All of the integrals in this expression are elementary; the important contribution comes 
from 




which diverges as 1/ y/T — T c near the Thouless curve. Since we have employed the (non self- 
consistent) T-matrix approximation, it is correct that this quantity diverges in Id; however, 
had we not handled this integration analytically (and instead proceeded with lattice sums 
for a large but finite lattice) we would have obtained a divergence proportional to 1/ (T — T c ) 
instead, which is incorrect in the thermodynamic limit. 
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The behaviour of the L — > oo (pseudo)gap A can now be identified through a (mean-field 
like) ansatz for the (pseudo)gap given by 

E(MO = . A \ • (29) 

This ansatz has been proposed by various authors; Moreo et a/.— remarked that the mean- 
field theory of Schrieffer et a/.— predicts such a form. The TPSC theory of Vilk and 
Tremblay^i leads to the same functional form for estimating the magnitude of the pseudo- 
gap. Also, the work of the group of Levin^ has produced and critiqued this ansatz. Lastly, 
we note that three of the present authors showed how such a form arises in a self-consistent 
T-matrix formulation.— 

We augment support for this ansatz by noting that when the (pseudo)gap parameter A 
is evaluated in the non self-consistent T-matrix approximation (in Id), focussing on the iQ n 
term in the Matusbara frequency sum, which diverges at T c , one expects that 

W v x°(g,Q) 

(5L V 1-IWM) 



A 2 - ^ E • OO) 



Here we note the following successes of this ansatz: (I) Using the finite-size scaling result of 
Eq. (|23|l, Eq. (jSH predicts that A 2 (L) = A 2 (L -> oo) + C(l/L 2 ), which also implies that 
a 1/L 2 scaling should be found for A; this is precisely the form that we find, as shown in 
Fig. |21 (II) If this ansatz is used to derive the DOS when the system is close to the Thouless 
criterion line for large L, good agreement is found. To be specific, if one solves for the 
temperature-dependent (pseudo)gap using Eq. (J30|) . and then uses this value in Eq. (J29|) . 
and then analytically continues to the real axis to find the spectral function A(k, u) through 
Eq. (JEj), one finds 



AiK ^ = K 1 + TfM ^ " v ^ i) + ^ " 7«M ^ + ' <31) 

Then, to facilitate a comparison with our partial fractions results, one must broaden the 
delta functions in Eq. (J3T)) in an analogous manner to that used in our partial fractions 
decomposition, that is each delta function is replaced by a Lorentzian with a fixed broadening 
of 0.12t. Finally, a direct comparison of the DOS from both the exact partial fractions 
expansion and the above described method yields very good agreement, in particular for the 
locations of the peaks closest to the "gapped region" . An example of such a comparison is 
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FIG. 4: A comparison of the "exact" DOS of a finite Id system, calculated within the NSC T- 
matrix approximation, to that predicted by Eq. I|29|). The (pseudo)gap parameter A is found from 
Eq. (|3U[). and then the density of states is calculated by using a wave vector sum of the spectral 
functions given in Eq. Q31|) . For these calculations the system parameters are L=128, U/t = —2, 
(n} = 0.5, and T/T C (L) = 1.0005, where a broadening of 0.12t has been used in both calculations. 

shown in Fig. |3J This agreement is exceptional given the crude nature of the approximation 
in Eqs. (EHUD- 

Thus, from Eq. (|25|l one expects that E diverges as 1/ (y/T — T c ), and then using Eq. (|2*9*Jl 
the (pseudo)gap should diverge as 1/(T — T c ) 1//4 (in fact, this result can be directly extracted 
from Eq. (|30)l). This is precisely the temperature dependence that we find - in Fig. El 
we show the L — ► oo extrapolated (pseudo)gap vs. the reduced temperature, denned by 
ST = (T — T c )/T c , fit to the form 1/5T K , and obtain an exponent k = 0.27 ± 0.03. 

Note that despite the potential difficulties of making a precise determination of the peak 
position (e.g., accounting for systematic errors introduced through the use of some fixed 
constant Lorentzian broadened peaks obtained from the partial fractions decomposition), 
the obtained fit to the data is excellent, with a critical exponent very close to the expected 
value of 1/4. 

We note that one may complete similar calculations in two dimensions, and again the 
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FIG. 5: The locations of the peaks, extrapolated to the thermodynamic limit using the previously 
described analysis, plotted as a function of the reduced temperature, defined by 5T = (T — T c )/T c , 
just above the superconducting transition temperature. The same parameters as in Fig.^are used. 
The solid line shows the least-squares fit of the data to the functional form A/t ~ (5T) K , where 
we find k = 0.27 ± 0.03. 

(pseudo)gap energy diverges to infinity as the transition temperature is approached from 
above. However, in two dimensions this divergence is logarithmic, a result that follows from 
Eq. ([HI) , and thus unless one is very close to the transition, this divergence is difficult to 
find. For example, the real-axis analysis of the T-matrix approximation in two dimensions 
of Metzner et a/.— finds spectral functions and pseudogaps similar to what we find, but 
they did not go to small enough reduced temperatures to observe the (unphysical) diverging 
pseudogap. (What is needed to suppress this unphysical divergence, which can be traced to 
the unphysical divergence of the effective vertex, is a theory beyond the non self-consistent 
T-matrix approximation . 31 ' 38 ) 

We believe that the sequence of calculations outlined above for one dimension, viz. (i) 
perform a partial fractions expansion of the pair propagator, thus obtaining the poles and 
residues of the one-particle Green's function (for a finite lattice) to arbitrary numerical 
accuracy, and (ii) extrapolating to the thermodynamic value of the (pseudo)gap as a function 
of the reduced temperature using Eq. (|23|). is justified, in part by the above demonstration. 
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IV. A PSEUDOGAP IN THREE DIMENSIONS 



We now describe the result that is central to the goal and message of this paper, viz. 
a single-particle DOS pseudogap above the superconducting T c can be present in three 
dimensions, even in the weak coupling limit. 

To begin, we first note the qualitative change that arises from an analysis of the behaviour 
of the self energy in three vs. one (or two) dimensions. That is, according to Eq. ([25)1 . even 
when the denominator of x{q.i 0) vanishes at q = the self energy remains finite. It is this 
difference that makes the study of the AHM on a three-dimensional lattice so attractive 
(no pun intended). That is, one correctly predicts a true phase transition at a nonzero 
temperature, and the Mermin- Wagner theorem does not disallow such a transition in three 
dimensions - it is the true superconducting phase transition that such a system will undergo. 

We have repeated the analysis of the previous section for a three-dimensional simple cubic 
lattice with periodic boundary conditions. As system parameters we have used \U\/t = 6 
and 3 (that is, \U\ = W/2 and \U\ = W/4, where W = 12t is the noninteracting bandwidth 
in three dimensions) . The first of these ratios is the same as the on-site attractive interaction 
energy to bandwidth that was used in the previous section's study of the Id AHM, while 
the second is in the weak-coupling limit. Also, in part due to the interest in the cuprate 
superconductors and doped Mott insulators, we have examined a (noninteracting) electronic 
density of (n)o = 1, viz. half filling. 

Results for the DOS for different size lattices are shown in Fig. El for \U\/t = 6, for a 
reduced temperature of 0.0001. As in the previous section, we have broadened the delta 
functions that emerge from the partial fractions decomposition with a fixed energy equal to 
3 % of the noninteracting bandwidth, which is 0.36t for the simple cubic lattice. It is to be 
emphasized that there are no true superconducting gaps in these spectra, in that the DOS 
is never reduced to zero at the Fermi level (which is located at u — 0). 

How does the DOS develop as a function of the attractive Hubbard energy U? In Fig.Elwe 
show the change of the DOS as \U\/t is decreased from 6 to 3, using a reduced temperature 
of 0.001 for each value of \U\, for a 12 3 lattice; that is, T/T C (\U\, L) = 1.001. As a reference 
curve, we also show the U = noninteracting DOS in the thermodynamic limit. 

The scaling of the peak positions, which we now label as A pg , vs. system size is shown 
in Fig. El and we again find a finite-size scaling behaviour that clearly goes as 1/L 2 . From 
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FIG. 6: The single-particle DOS for simple cubic lattices of size 6 3 , 8 3 , 10 3 , and 12 3 , for \U\/t = 6 
and 1/2 filling, at a temperature of T/T C (L) = 1.0001. We have employed a fixed broadening of 
3% of the noninteracting bandwidth, namely 0.36t. 

this figure one sees that all extrapolated values for small reduced temperatures are very 
close to one another (approximately Q.5t for these material parameters), and do not go to 
zero, or infinity, as the reduced temperature is lowered towards zero. Further, one can use 
BCS theory^ to calculate the superconducting energy gap at T = 0, and for 1/2 filling and 
\U\/t = 6 one finds a gap of about 2.3t, roughly 1/3 of the energy of the pseudogap found 
in Fig. |S1 (which is evaluated just above T c , and not T = 0!). 

How does the pseudogap behave as a function of temperature? As seen in Fig. |HJ for 
5T = 0.1 the pseudogap is roughly the same as its value at lower temperatures, in the ther- 
modynamic limit. In fact, for a temperature of 1.5T C we find a pseudogap that is about 85% 
of its value as ST — ► 0. So, the pseudogap energy does not change appreciably with temper- 
ature for small reduced temperatures (less than 0.1). At sufficiently high temperatures the 
pseudogap disappears, but we find that it does so by the DOS filling in, as opposed to the 
pseudogap energy going to zero. 

As to whether or not the pseudogap persists for lower couplings, we have completed a 
similar analysis for \U\/t = 3 at half filling, and again find evidence for a pseudogap - the 
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FIG. 7: DOS for \U\/t = 6 and 3 at half filling, for a 12 3 simple cubic lattice, at a temperature of 
T = 1.001 T c . In this plot we have used the same broadening (0.36t) for both [7^0 curves. As 
a reference curve we also show the DOS of the noninteracting system in the thermodynamic limit. 

analogue of Fig. |H1 is shown in Fig. We find the striking result that the extrapolated 
gap for smaller reduced temperatures is within 5% of the BCS gap! We emphasize that this 
matching of energy "gaps" is found between a T = BCS gap and a T c < < 1.001T C 
pseudogap. 

Our result, that is the existence of a pseudogap in a model of a three-dimensional strongly 
correlated electronic system, agrees with the earlier (three-dimensional) work of Levin et 
al.M However, they propose that the pseudogap they find can be understood in terms of a 
T-matrix resonance caused by intermediate to large U, whereas we are finding the full value 
of the BCS gap above T c even for weak coupling. 

For smaller reduced temperatures and large lattices our partial fractions algorithm fails, 
due to the high accuracy that must be enforced to find converged pole locations and residues. 
To be specific, we were unable to find a converged partial fractions decomposition even after 
30 days of computation for L — 12 for reduced temperatures less than 0.0005. 

Clearly, a thorough search of parameter space, that is for other electronic densities and 
for other ratios of \U\/t, is called for. Unfortunately, as mentioned earlier, it takes roughly 10 
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FIG. 8: The peak locations of the DOS, labelled as A pg and expressed in units of t, for 3d simple 
cubic lattices of size 6 3 , 8 3 , 10 3 , and 12 3 , for \U\/t = 6 and 1/2 filling, plotted as a function of one 
over the the system's linear size (Li) squared, for a variety of temperatures (shown in the legend) 
just above the superconducting transition temperature. The solid lines are least squares fits for 
each temperature. The T = superconducting gap (in the thermodynamic limit), found from the 
solution of the BCS gap equation, is also shown. 

days to get each data point for a 12 3 lattice (for one temperature and one chemical potential), 
so new implementations of this approach are being explored in the hope of finding a more 
expedient algorithm. 

V. CONCLUSIONS: 

We have introduced a technique that allows for the extraction of the real-axis dynamical 
properties of the pair susceptibility, the self energy, and the spectral functions, for a strongly 
correlated electronic system when unrenormalized propagators are used. We have applied 
this technique to the attractive Hubbard model, and examined the properties of the density 
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FIG. 9: The same quantities as in Fig. El are plotted for \U\/t = 3. Note that the extrapolated 
values of the pseudogap, at temperatures just above the superconducting transition temperature, 
are within 5% of the T = BCS gap. 

of states in the thermodynamic limit. Of particular interest for us is three dimensions, since 
a true transition at a nonzero temperature to a superconducting phase should be present 
for this model of correlated electrons, and for weak coupling, described quite well by the 
non self-consistent T-matrix approximation. We find a pseudogap in three dimensions that 
persists in the weak coupling limit, and whose energy is comparable to the T = BCS gap 
energy in the same weak coupling limit; for larger \U\/t we find a pseudogap energy that is 
much larger than the T = BCS gap energy . 

Acknowledgments 

We wish to thank Andrew Callan- Jones and Andre-Marie Tremblay for useful discussions. 
This work was supported, in part, by the NSERC of Canada. One of us (FM) acknowledges 
the support of the Quantum Materials Program of the CIAR, and ICORE (Informatics 
Circle of Research Excellence) of Alberta. 



25 



1 J. Bardeen, L.N. Cooper and J.R. Schrieffer, Phys. Rev. 106, 162 (1957); Phys. Rev. 108, 1175 
(1957). 

2 P.W. Anderson, Phys. Rev. 112, 1900 (1958). 

3 N.N. Bogoliubov, N.V. Tolmachev, and D.V. Shirkov, A New Method in the Theory of Super- 
conductivity, Consultants Bureau, Inc., New York (1959). 

4 G. Rickayzen, Phys. Rev. 115, 795 (1959). 

5 A. Bardasis and J.R. Schrieffer, Phys. Rev. 121, 1050 (1961). 

6 A.B. Migdal, Zh. Eksp. Teor. Fiz. 34, 1438 (1958) [Sov. Phys. JETP 7, 996 (1958)]. 

7 CM. Eliashberg, Zh. Eksperim. i Teor. Fiz. 38 966 (1960) [Soviet Phys. JETP 11 696 (I960)]. 

8 See the various articles in Superconductivity, edited by R.D. Parks (Marcel Dekker, New York 
(1969)). 

9 V.L. Ginzburg and L.D. Landau, Zh. Eksp. Teor. Fiz. 20, 1064 (1950). 

A.S. Alexandrov and N. Mott, Polarons and Bipolarons (World Scientific, Singapore, 1995), 
and references therein. 

1 D. Feinberg, S. Ciuchi, and F. de Pasquale, Int. J. Mod. Phys. B 4 1317 (1990). 

2 D. Emin, Phys. Rev. B 48, 13691 (1993). 

3 J. Ranninger, in Bose-Einstein Condensation (Cambridge University Press, 1995), edited by A. 
Griffin, D.W. Snoke, and S. Stringari, p. 393. 

4 V.J. Emery and S.A. Kivelson, Nature 374, 434 (1995). 

5 M. Randeria, in Bose-Einstein Condensation (Cambridge University Press, 1995), edited by A. 
Griffin, D.W. Snoke, and S. Stringari, p. 355. 

6 P.B. Allen, in Modern Trends in the Theory of Condensed Matter, Proceedings, Karpacz, 1979 
(Springer Verlag, New York, 1980), edited by A. Pekalski and J. Przystawa, p. 388. 

7 V. Ambegaokar, in Superconductivity, edited by R.D. Parks (Marcel Dekker, New York (1969)) 
Vol. 1, p.259. 

8 D.J. Thouless, Ann. Phys. 10, 553 (1960). 

9 K. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 

:0 A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems (McGraw-Hill, 
1971). 



26 



21 L.P. Kadanoff and P.C. Martin, Phys. Rev 124, 670 (1961). 

22 G. Baym and L.P. Kadanoff, Phys. Rev. 124, 287 (1961). 

23 G. Baym, Phys. Rev. 127 1391 (1962). 

24 A. Schmid, Z. Physik 231, 324 (1970). 

25 S. Marcelja, Phys. Rev. B 1, 2351 (1970). 

26 B.R. Patton, Phys. Rev. Lett. 27, 1273 (1971); PhD thesis (Cornell), 1970. 

27 V.M. Loktev, R.M. Quick, and S.G. Sharapov, Phys. Repts. 349, 2 (2001). 

28 B. Janko, J. Maly, and K. Levin, Phys. Rev. B 56, 11407 (1997); I. Kosztin, Q. Chen, B. Janko, 
and K. Levin, Phys. Rev. B 58, 5936 (1998). See also Q. Chen, K. Levin, and I. Kosztin, Phys. 
Rev. B 63, 184519 (2001), and references therein. 

29 N.E. Bickers, D.J. Scalapino, and S.R. White, Phys. Rev. Lett. 62, 961 (1989); N.E. Bickers 
and D.J. Scalapino, Ann. Phys. (N.Y.) 193, 206 (1989); N.E. Bickers and S.R. White, Phys. 
Rev. B 43, 8044 (1991). 

30 S. Verga, F. Marsiglio, and R.J. Gooding (to be submitted). 

31 J. Vilk and A.-M.S. Tremblay, J. de Physique I (France) 7, 1309 (1997). 

32 A.J. Leggett, J. de Physique, C7, 41, 19 (1980); 

A.J. Leggett, in Modern Trends in the Theory of Condensed Matter, edited by S. Pekalski and 
J. Przystawa (Springer, Berlin, 1980) p. 13. 

33 P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 195 (1985). 

34 S. Schmitt-Rink, C. Varma, and A.E. Ruckenstein, Phys. Rev. Lett. 63, 445 (1989). 

35 J.M. Luttinger, Phys. Rev. 119, 1153 (1960). 

36 J. Serene, Phys. Rev. B 40, 10873 (1989). 

37 A. Tokumitu, K. Miyake, and K. Yamada, Phys. Rev. B 47, 11988 (1993). 

38 K.S.D. Beach, R.J. Gooding, and F. Marsiglio, Phys. Lett. A 282, 319 (2001). 

39 R. Fresard, B. Glaser and P. Woffle, J. Phys. Cond. Mat. 4, 8565 (1992). 

40 R. Haussmann, Z. Phys. B 91, 291 (1993); R. Haussmann, Phys. Rev. B 49, 12975 (1994). 

41 M. Y. Kagan, R. Fresard, M. Capezzali, and H. Beck, Phys. Rev. B 57, 5995 (1998). 

42 B. Kyung, E. G. Klepfish, and P. E. Kornilovitch, Phys. Rev. Lett. 80, 3109 (1998). 

43 V.P. Gusynin, V.M. Loktev, and S.G. Sharapov, JETP 88, 685 (1999). 

44 J.R. Engelbrecht and A. Nazarenko, Europhys. Lett. 51, 96 (2000). 

45 D. Rohe and W. Metzner, Phys. Rev. B 63, 224509 (2001). 

27 



46 Y. Yanase and K. Yamada, J. Phys. Soc. Jpn. 70, 1659 (2001). 

47 J. Quintanilla, B.L. Gyorffy, J.F. Annett, and J. P. Wallington, Phys. Rev. B 66, 214526 (2002). 

48 P. Pieri an G.C. Strinati, Phys. Rev. B 61, 15370 (2000). 

49 S. Allen and A.-M.S. Tremblay, Phys. Rev. B 64, 075115 (2001). 

50 B. Kyung, S. Allen, and A.-M.S. Tremblay Phys. Rev. B64, 075116 (2001). 

51 We make this statement with the understanding that we are not attempting to discuss 
Kosterlitz-Thouless physics in this paper. 

52 T. Timusk and B.W. Statt, Rep. Prog. Phys. 62, 61 (1999). 

53 G. Rickayzen, Green's Functions and Condensed Matter (Academic Press, Toronto, 1980). 

54 H. Fukuyama, O. Narikiyo, and Y. Hasegawa, J. Phys. Soc. Jpn. 60, 372 (1991); H. Fukuyama, 
Y. Hasegawa, and O. Narikiyo, J. Phys. Soc. Jpn. 60, 2013 (1991). 

55 F. Marsiglio and R.J. Gooding (unpublished). 

56 This terminology was first used in T. Gherghetta and Y. Nambu, Phys. Rev. B 49, 740 (1994). 

57 Due to several important technicalities, e.g., to ensure that we purge spurious double poles 
from any expansion generated by this algorithm, a very high numerical accuracy is required. 
We have coded in Maple Vr5.1 to achieve the extraction of these numbers, and have used a 
relative accuracy as high as 10 -150 , noting that such a number can be required for converged 
and accurate pole and residue evaluations. 

58 N.D. Mermin and Wagner, Phys. Rev. Lett. 17, 1133, 1307 (1966); PC. Hohenberg, Phys. Rev. 
158, 383 (1967). 

59 A. Moreo, D.J. Scalapino, R.L. Sugar, S.R. White, and N.E. Bickers, Phys. Rev. B 41, 2313 
(1990). 

60 J.R. Schrieffer, X.-G. Wen, and S.C. Zhang, Phys. Rev. B 39, 11663 (1989). 

61 We note that to complete the \U\/t = 3 analysis we have had to use a smaller broadening 
parameter of 0.12t, in contrast to the value of 0.36t used for \U\/t = 6. This is because the 
energy scales are much smaller, and a broadening of 0.36i is in fact more than 50% of the BCS 
gap; to properly resolve the energies for smaller \U\/t we thus used a smaller broadening. 



28 



